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ABSTRACT 

Using a sample of nearly 700 quasars with strong (Wo (2796) > 0.8 A) Mgll absorption lines detected 
in the Early Data Release of the SDSS, we demonstrate the feasibility of measuring the photometric 
properties of the absorber systems by stacking SDSS imaging data. As Mgll lines can be observed in 
the range 0.37 < z a b s < 2.2, the absorbing galaxies are in general not identified in SDSS images, but 
they produce systematic light excesses around QSOs which can be detected with a statistical analysis. 
In this Letter we present a 6 a detection of this effect over the whole sample in i-band, rising to 9.4 a 
for a low-redshift subsample with 0.37 < z a b s < 0.82. We use a control sample of QSOs without 
strong Mgll absorption lines to quantify and remove systematica with typical 10-20% accuracy. The 
signal varies as expected as a function of absorber redshift. For the low subsample we can reliably 
estimate the average luminosities per Mgll absorber system in the g, r, and i bands and find them to 
be compatible with a few-hundred-Myr old stellar population of M r ~ —21 in the rest frame. Colors 
are also consistent with typical absorbing galaxies resembling local Sb-c spirals. Our technique does 
not require any spectroscopic follow-up and does not suffer from confusion with other galaxies arising 
along the line-of-sight. It will be applied to larger samples and other line species in upcoming studies. 
Subject headings: quasars: absorption lines - galaxies: photometry - galaxies: statistics - techniques: 
photometric 



INTRODUCTION 

A number of studies have inferred and then shown 
that strong metal absorption lines in a quasar spec- 
trum are due to gaseous clou ds associated to galaxies 
close to the line-of-sight (e.g. iBahcall fc Spitzerlll969t 
[ Bergeron fc Boissel Il991t ISteidel. Dickinson, fc Perssonl 
1994). However, the systematic characterization of the 
emission properties of these absorber systems is still 
limited, especially a t high redshifts (see for example 
iChurchill et"aHl200l . In order to use absorption-line 
measurements for investigating galaxy formation and 
evolution, characterizing the underlying population of 
galaxies seen in absorption is an important task to 
achieve. 

Among metal absorption lines, the Mgll doublet 
AA2796, 2803, detectable from the ground at 0.2 < z < 
2.2, has been widely used. It arises in gas spanning more 
than five decades of neut ral hydrogen column density 
(e.g. IChurchill et alJfeOOOj) and the strong systems, de- 
fined by W (AA2796) > 0.3 A are found in the vicin- 
ity of galaxies spanning a large range of morphologies. 
Previous attempts to constrain the properties of Mgll 
absorbers have relied on deep imaging and spectroscopic 
follow-up of (arbitrarily faint) galaxies in the quasar field 
in order to identify a potential galaxy responsible for the 
absorption seen in the quasar spectrum. Such expensive 
studies have been limited to samples of a few tens of 
quasars. 
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Here we propose a new approach to measure the sys- 
tematic photometric properties of large samples (N » 
100) of absorbing systems, which combines statistical 
analysis of both spectroscopic and imaging datasets of 
the S loan Digital Survey (hereafter SDSS, lYork et all 
2000;. As Mgll lines can only be detected in SDSS spec- 
tra in the range 0.37 < z < 2.2, the absorbing galax- 
ies are not identified and, in general, only marginally 
detected in individual SDSS images. However, they do 
produce systematic light excesses around the background 
QSOs that can be detected with a statistical analy- 
sis consisting of stacking a large number of absorbed 
QSO images. Such a statistical approach has been al- 
ready proved to be very powerful in detecting and char- 
acterizing the v ery low-surface brig htness diffuse light 
in galaxy halos |Zib etti et al.ll2004|) and the intraclus- 
ter light (jZibetti et alJ l2005l. In this Letter we demon- 
strate that significant detections can be obtained with 
this technique also for localized sources (the absorbing 
galaxies) in the presence of strong background contami- 
nation. Preliminary results on the photometric proper- 
ties of the low z (< 0.82) absorbers are also analyzed. 
More detailed results and analysis for a larger sample 
will be presented in a forthcoming paper. 

The paper is organized as follows. In Section ^ wc 
present the samples of QSOs used in our analysis. The 
image processing and stacking are described in Sec. 
The results are presented and analyzed in Sec. Eland con- 
clusive remarks and forthcoming developments are out- 
lined in Sec. 01 Standard ACDM cosmology (fi m = 0.3, 
Oa = 0.7) and H — 70 km sec -1 Mpc -1 are adopted 
throughout the paper. 



1. THE DATA 



1.1. 



The Mgll absorber sample 

We use th e sample of Mgll a bsorption lir 
compiled by iNestor et al.l l)2004|) and based 
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Fig. 1. — Redshift distribution of Mgll absorbers with Wq > 
0.8 A (shaded histogram), of their background quasars (solid line), 
and of the reference quasars (dotted line). 



EDR data <|Stoughton et al.ll2002fl . We focus our analy- 
sis only on strong systems with a rest equivalent width 
Wo(2796) > 0.8 A. In this range the detection complete- 
ness is above 95%. In this section we briefly summarize 
the main steps involved in the absorption line detection 
procedure. For details on t he quasar and absorber cat- 
alogs, we refer the reader to lSchneider et all l|2002fl and 
iNestor et all (|2004|) . 

The SDSS EDR provides approximately 3700 QSO spec- 
tra of sufficiently high redshift to allow the detection 
of intervening Mgll absorption lines. The continuum- 
normalized spectra were searched for Mgll AA2796, 2803 
doublets using automated routines and interactive confir- 
mation/rejection of candidate systems. The rest equiva- 
lent widths Wq of confirmed systems were measured us- 
ing an optimal extraction method employing Gaussian 
line-profiles (or more complex ones where appropriate). 
All systems were checked for blends and other special 
cases. A 5a significance level was required for all A2796 
lines, as well as a 3a significance level for the corre- 
sponding A2803 lines. Only systems 0.1c blue- ward of 
the quasar redshift and red-ward of Ly— a emission (so 
as to avoid the Lyman-a forest) were accepted. Systems 
with separations of less than 500 km/s were considered 
as a single system. In order to avoid confusion originat- 
ing from the superposition of multiple absorbers along 
the same line-of-sight, we exclude from our sample all 
QSOs where distinct W (2796) > 0.8 A absorber sys- 
tems are detected: they represent 16.3% of QSOs with 
Wo (2796) > 0.8 A absorptions. The final sample com- 
prises 683 QSOs, out of which 2 are excluded because 
they lie too close to bright stars. The absorber and QSO 
redshift distribution are shown in Figure ^ 

1.2. The reference quasar sample 

In order to measure possible systematic light excess 
around QSOs that could contaminate the signal origi- 
nating from absorbing galaxies, we use a control sam- 
ple made of QSOs without strong absorbers. For each 
QSO with an absorber system, we randomly look for 
another QSO without any detected absorber down to 
Wq 2796 = 0.8 A, such that their redshifts are matched 
within 0.1 5 and i-band brightnesses differ by less than 
0.5 mag. In addition, we require that their spectra have 
similar signal-to- noise (S/N) at the wavelengths corre- 
sponding to the Mgll absorption features. The last two 
conditions also imply that the average overall S/N per 

5 Only in one case the best possible match is within 0.3. 



pixel in the spectra of the two samples differs by less than 
0.1. Thus our procedure ensures that (i) the probability 
of contamination by undetected Mgll systems along the 
line-of-sight is the same for both QSOs and (ii) photomet- 
ric systematics which may possibly arise from residuals 
of the PSF subtraction (see Sec. I2.1l andl3"|l affect both 
samples in the same way. Point (i) has also been ver- 
ified by means of Monte Carlo simulations. The QSOs 
selected in this way are called reference QSOs in the fol- 
lowing. Their redshift distribution is also shown in Fig. 
^ (dotted histogram). 

2. IMAGE STACKING 

The image stacking technique aims at producing a 
high S/N average image of the QSOs, possibly sur- 
rounded by absorbing galaxies. It consists of four ba- 
sic steps: (i) shifting to make images superposable; (ii) 
intensity rescaling to a uniform photometric calibra- 
tion; (iii) masking of all unwanted sources; (iv) average- 
combination of images, after rejecting masked regions. 
All the procedures will be thoroughly described in a 
forthcoming paper. Here we summarize the main points. 

We use the raw "corrected frames" (fpC, bias and flat- 
field corrected) pub licly available in the T hird Data Re- 
lease of the SDSS ijAbazaiian et alJ [2005. DR3). Steps 
(i) and (ii) rely on the centroid calculations and on the 
photometric calibrations provided by the DR3 catalogs 
and are performed using standard IRAF tasks. A bi- 
cubic spline interpolant is used to shift images. Regard- 
ing step (iii), we note that an ideal mask should retain 
only the QSO and the absorbing galaxy in each image, 
while blanking out completely all other sources in order 
to minimize the background noise. In practice we de- 
sign our masking algorithm to leave the QSO unmasked 
and to blank all sources which are likely to be in the 
foreground with respect to the absorber. For each QSO 
we create three independent masks in g- 1 r- and i-band, 
that are eventually combined in the final mask. In each 
band, a first mask is obtai ned from the segmentatio n im- 
age output by SExtractor ijBertin k, Arnoutslll996|) for a 
clean image of the field, from which we have subtracted 
the de-blended "Atl as Image" of the Q SO produced by 
the PHOTO pipeline ijLupton et al.ll2001|) . This step en- 
sures optimal de-blending of sources from the QSO and 
leaves the QSO unmasked. We retain the masks for 
all objects with a flux exceeding that expected from 
iBruzual &: Chariot s population synthesis models (|2003L 
hereafter BC03) for an extreme case of a luminous blue 
galaxy: a stellar population with M< = —22.0 produced 
by a burst 10 7 Gyr long, observed immediately after the 
burst end at the redshift of the absorber. Fainter objects 
are left unmasked. Our tests have shown that the final 
results do not critically depend on the assumed spec- 
tral energy distribution. Finally, the masks in the three 
bands are combined and conservative masks are applied 
to all objects classified as stars by PHOTO with r < 20.0. 

After subtracting the background value evaluated on 
the masked image in an annulus between 200 and 250 pix- 
els, the original images are shifted, intensity-calibrated, 
and averaged, excluding all masked pixels. Azimuthally 
averaged surface brightness (SB) profiles are extracted 
in a series of circular apertures which are optimized to 
integrate on larger areas at larger radii. A more accu- 
rate background level for these profiles is recomputed 
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on the stacked image between 100 and 150 pixels. The 
complete covariance matrix of the extracted background- 
subtracted SB profile is evaluated using the jackknife 
technique. 

2.1. Subtracting the effective PSF 

To isolate the excess light due to intervening absorbers, 
the mean point spread function (PSF) must be sub- 
tracted. The PSF depends on several parameters, such 
as the time of the observation, the position on a given 
camera column, and the color of the object. In order 
to build an effective PSF which is representative of the 
flux-averaged PSF of our QSOs extended to large an- 
gular distances, we stack stars selected to match each 
individual QSO. We first require stars to be observed in 
the same observing run and same camera column as the 
QSO. These stars are then ranked to minimize a com- 
bination of the following quantities: (i) the difference of 
Gaussian FWHM with respect to the QSO; (ii) the an- 
gular separation from the QSO; (iii) the differences in 
colors with respect to those of the QSO. To map the 
PSF with reasonable S/N at large distances, bright un- 
saturated stars are also preferred (r < 17.0 mag is always 
required). For each QSO in each band, the star that best 
fulfills all these criteria is chosen. To obtain the correct 
flux-weighted average PSF, the intensity of each star is 
rescaled to the same intensity of the corresponding QSO 
before stacking. 

3. RESULTS 

In Fig. we show the SB residuals measured in the 
i-band after the PSF subtraction for the sample of 681 
QSOs with absorbers (upper panel) and for the refer- 
ence QSO sample (lower panel). The stacked PSFs are 
normalized with respect to the QSO profiles within an 
aperture of 4 pixels (1-6"). The vertical scale is linear 
in flux intensity and the horizontal solid line is the back- 
ground level. The left vertical axis displays the surface 
brightness in mag arcsec~ 2 (in excess or deficit). Er- 
ror bars are computed from the jackknife covariance ma- 
trix and include the budget from the background uncer- 
tainty. Central normalization uncertainties are found to 
be negligible. We note that within ~ 4" both profiles 
display SB excesses: these systematic departures from 
the stellar PSF is due to reduced efficiency of source de- 
blending close to the QSO and, partly, to the presence of 
the QSO host galaxy. However, we see very clearly that 
from 3" outwards the SB excess is significantly larger 
for absorbed QSOs. In particular, between 4 and 8" the 
absorbed QSOs display positive SB excesses that are sig- 
nificant at 2 to 4 a level each, while the profile of the 
reference QSOs is fully consistent with the background 
level. Similar results are found for the r- and g-band 
as well. Note that galaxies which produce undetected 
absorbing systems at different rcdshifts are expected to 
be present with the same statistically uniform surface 
density as in the background. Therefore the SB excess 
in the absorbed with respect to the reference QSOs are 
solely due to galaxies related to the selected absorption 
systems. 6 

6 We cannot exclude, however, that more than one galaxy per 
absorber contributes to the signal. This might be the case if ab- 
sorbing galaxies preferentially reside in groups. 
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Fig. 2. — i-band QSO— PSF residuals as a function of radius R 
for the absorbed QSOs (upper panel) and the reference unabsorbed 
QSOs (lower panel). The y-axis scale is linear in flux intensity, as 
indicated on the right-side (arbitrary count units per pixel). On 
the left the corresponding levels of the surface brightness excess 
or deficit are reported in mag arcsec - 2 . The horizontal solid line 
indicates the background level. Error bars are computed with a 
jackknife algorithm. 



In Table we report the fluxes obtained by integrat- 
ing the PSF-subtr acted SB within 1.6 and 10" (~ 80 kpc 
in the z a b s interval considered) for the absorbed QSOs 
in g-, r-, and i-band (columns 1, 3, and 5, respectively, 
in mag). The attached random errors (1 a) are com- 
puted by propagating the errors on the residual SB pro- 
files using the complete jackknife covariance matrix. In 
columns 2, 4, and 6 we report the systematic shifts that 
must be applied to the raw magnitudes integrated on 
the absorbers samples to correct for the bias in the in- 
ner 5", as measured from the reference sample. Corre- 
sponding 1 a upper and lower limits for this correction 
are also attached. Measurements are reported for the 
whole sample (681 QSOs) and for three subsamples of 
absorber systems binned in z a ^ s : 0.37 < z a b s < 0.82, 
0.82 < z a b s < 1.26 and 1.26 < z a b s < 2.3. Assuming the 
best estimated value for the systematic corrections, for 
the complete sample we have ~ 4er detection in g- and 
r-band, and 6cr in i. However, even with the most con- 
servative corrections our detections are significant at 2.9 
(g and r) and 5.3 a (i). As expected, looking at the three 
redshift bins, we note that the signal is clearly dominated 
by the low-z a b s one. For this subsample the significance 
is 2.6 a in g, 5.6 in r, and 9.4 in i. In the two higher 
redshift bins the detections in individual bands are not 
statistically significant, however the best estimates we 
derive are consistent with the dimming determined by 
the larger distance. 

We now focus on the low-z a b s subsample to derive more 
physical information about the absorbing galaxies. In 
order to check whether the signal is dominated by a few 
anomalously bright objects, we have computed the lu- 
minosity function (LF) of all resolved sources below the 
masking threshold around the absorbed QSOs and sta- 
tistically subtracted the analogous LF in the reference 
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TABLE 1 

Integrated magnitudes of QSO absorbers (see Section 0. 



Sample 


g 


r 


i 




m(< 10") Am sys. 


m(< 10") Am sys. 


m(< 10") Am sys. 




(1) (2) 


(3) (4) 


(5) (6) 



All 23.90;«;« +0.13+g;JJ 22.98-° it +0-^+° T 5 22.31;°;^ +0.287.g : °8 

* abs <0.82 23.76;0;2| +0.00^°;°° 22.257°;^ +0.10+";!° 21.407_°;°f +0.157.°;°5 

0.82 < 2abs < 1.26 24.55;°;« +0.17;°;^ 23.517.^ +1.03+%£ 23.63;°;" +0.74;?;^ 

* abB >1.26 23.63+g;gg +0A5-° 3 7 l 23.69^°;^ +1.21-°;g| 23.09;gJ« +0727^ 



QSOs fields. We find that less than 5% of the flux is con- 
tributed by galaxies that are more than 2 mag brighter 
than the average values given in Table while 80% is 
contributed by objects brighter than the average. The 
properties we derive are thus well representative of the 
luminosity weighted average galaxy emission associated 
to absorbing systems. 

From the photometric data presented above we derive 
colors within 10" aperture: g - r = 1.41 ± 0.35±g;if and 
r — i = 0.80 ± O.IG^q'^, (the first uncertainty is ran- 
dom error on the direct measure, the second accounts 
for errors on the systematic correction). By neglecting 
the distribution of absorbers in redshift and assuming 
the median redshift of the sample, the comparison with 
BC03 models of solar-metallicity dust-free single stellar 
populations indicates that the luminosity-weighted av- 
erage SED of Mgll absorbing galaxies is dominated by 
fairly young stellar populations (~ 500 Myr old), with 
a rest-frame r-band absolute magnitude of ~ —21. The 
inferred rest-frame B — K value is about 2.5 mag, which 
is compa tible with although slightly bluer than the value 
found bv lSteidel. Dickinson, k, Perssonl |l994) using di- 
rect imaging and spectroscopic follow up. A more sophis- 
ticated analysis, in which we compare the observed colors 
with the z-folding o f UV-optical galaxy templates (from 
iKinnev et aljfl99q) normalized to a fixed Mi, confirms 
that the "average" absorbing galaxy has SED similar to 
local intermediate spirals (Sb-c). More details will be 
presented in a forthcoming paper. 

4. CONCLUSIONS AND OUTLOOK 

In this Letter we have used a sample of 681 SDSS- 
EDR QSOs with Mgll absorption lines (W > 0.8 A) 
to demonstrate the feasibility of studying the statisti- 
cal photometric properties of galaxies associated to Mgll 
absorbing clouds with the SDSS. By stacking images of 
quasars and subtracting the PSF, we have obtained sig- 
nificant detections (up to 9.4 a in i-band for the low- 
redshift subsample) of surface brightness excesses around 
QSOs with strong Mgll absorbers in the g, r and i bands. 
Systematic biases are measured and removed with 10 to 
~ 20% accuracy (for the total and low-,z a bs samples) by 
using a control sample of QSOs without strong absorbers. 
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As expected, the measured excess light decreases as 
a function of absorber redshift. At z ~ 0.6, the mean 
colors {g — r = 1.4 and r — i = 0.8) indicate the pres- 
ence of relatively young stellar populations and are con- 
sistent with the typical absorbing galaxy resembling a 
local Sb-c spiral. The average luminosity per absorber 
within ~ 80 kpc impact parameter is similar to L*. 
Both results are compatible with previous finding by 
iSteide l. Dickin son, fc Perssonl 1^)94). A detailed analy- 
sis of impact parameter distributions and constraints on 
the corresponding stellar populations will be presented 
in a forthcoming paper. 

Previously, the largest photometric study of Mgll ab- 
sorbing galaxies was limited to ~ 60 low redshift systems 
ijSteidel. Dickinson. Ez Perssonlll994j) due to the expen- 
sive deep imaging and spectroscopic follow up. Our tech- 
nique allows the study of much larger samples: in this 
Letter we have analyzed a sample more than ten times 
larger, and it will be soon applied to the SDSS absorber 
database (York et al., in preparation), i.e. another order 
of magnitude larger. Although the statistical nature of 
our method makes it very difficult to study the real in- 
trinsic variance of the absorbing galaxy properties, it can 
provide robust statistical results and also minimize the 
problem of confusion due to chance line-of-sight superpo- 
sition of galaxies. Moreover this technique can be applied 
to different absorption line species (Fel, Fell, CIV, DLAs, 
etc.) to measure their impact parameters and investigate 
their spatial distribution in the gas around galaxies. Fi- 
nally, it can also be used to investigate the properties of 
intrinsic absorber systems and quasar host galaxies. 
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